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In this paper we investigate the limits of control for mixed-state quantum systems. The constraint 
of unitary evolution for non-dissipative quantum systems imposes kinematical bounds on the opti- 
mization of arbitrary observables. We summarize our previous results on kinematical bounds and 
show that these bounds are dynamically realizable for completely controllable systems. Moreover, 
we establish improved bounds for certain partially controllable systems. Finally, the question of 
dynamical realizability of the bounds for arbitary partially controllable systems is shown to depend 
on the accessible sets of the associated control system on the unitary group U(N) and the results 
of a few control computations are discussed briefly. 



I. INTRODUCTION 

Recent advances in laser technology have opened up 
new possibilities for laser control of quantum phenomena 
such as control of molecular quantum states, chemical re- 
action dynamics or quantum computers to mention only 
a few. The limited success of initially advocated con- 
trol schemes based largely on physical intuition in both 
theory and experiment [jl), has prompted researchers in 
recent years to study these systems systematically using 
control theory. One of the many interesting questions 
that has arisen is the issue of dynamical realizability of 
bounds on the optimization of observables imposed by 
kinematical contraints, which we shall address in this pa- 
per. 

The answer to this question is not just of theoretical in- 
terest, it is a matter of practical importance as well since 
most algorithms designed to find optimal controls ^ |j| 
to drive a system are based on a set of differential equa- 
tions providing necessary but not sufficient conditions for 
a global maximum, i.e., these algorithms may produce 
controls that steer the system to a local maximum or 
minimum, for instance. Independent knowledge of dy- 
namically attainable bounds on the expectation value of 
the observable makes it possible to determine the effec- 
tiveness of a given control in achieving the control objec- 
tive of maximizing or minimizing the expectation value 
of a given observable. 



II. QUANTUM STATISTICAL MECHANICS 
MODEL 



We consider a quantum-mechanical system whose pure 
states form a separable Hilbert space H. A mixed state 
is an ensemble of orthonormal pure quantum states 
with a discrete probability distribution assigning each 
pure state a certain probability Wk such that < Wk < 1 
and J2k w k — 1- Such a state can always be represented 



by a density operator p on TL with eigenvalue decompo- 
sition 



(1) 



where Wk are the eigenvalues, and the corresponding 
normalized eigenstates of p. Unless otherwise specified, 
we shall use the word "state" in the following to refer to 
a mixed quantum state represented by a density operator 
P- 

If the system is Hamiltonian, the time-evolution of p 
is given by 

p(t) = U{t,t )poU(t,t Q )\ (2) 

where po = p(to) and C/(i,i ) is the time-evolution 
operator of the system satisfying the time-dependent 
Schrodinger equation 

d 



ihj-U(t,t ) 
at 



HU(t,to). 



(3) 



H is the total Hamiltonian of the system. 

Observables are represented by Hermitian operators A 
on TL and we define their expectation value to be the 
ensemble average 



(A(t)) = Tr (Ap(t) 



(4) 



The aim of controlling the system is to maximize the 
expectation value of a chosen observable A at a certain 
target time tp by driving the system using a set of opti- 
mal control fields f(i) = (fi(t), ■ ■ ■ , /m(*))- If the control 
fields are sufficiently weak, it can be assumed that the 
system is control-linear, 



H = H 



M 

E 



(5) 



where Hq is the internal Hamiltonian of the system and 
H m , m > 1, is the interaction Hamiltonian for the field 

fm- 
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III. KINEMATIC AL RESTRICTIONS ON THE 
DYNAMICS 

If dissipative effects are negligeable then the quan- 
tum system is Hamiltonian and therefore the evolution of 
the system has to be unitary no matter how the system 
is driven. This observation has profound implications. 
Given an initial mixed state po, the only kinematically 
attainable target states p(t F ) are 

p{t F ) = U(t F ,t )poU(t F ,to)\ (6) 

where U{t Fl to) is a unitary transformation. 

This kinematical constraint on the dynamical evolution 
leads to bounds on the expectation value of arbitrary ob- 
servables that arc completely independent of the control 
functions and thus impose kinematical restrictions on the 
optimization of observables, summarized by the following 
two theorems 

Theorem 1 Let A be a Hermitian operator on 7i with 
eigenvalue decomposition 

m 

A = J2aJ( ai ), (7) 

i=l 

where /(flj) is the projector onto the eigenspace E(ai), 
and let Xi > X2 > ■ ■ ■ > Xn be the eigenvalues a.i, counted 
with multiplicity, and ordered in a decreasing sequence. 
Then we have 

N N 

^ X N - k +iWk < Tr (Ap(t)j < ^ ^kw k , (8) 

fc=l k=l 

provided that the weights Wk are ordered in a decreasing 
sequence, i.e., ui\ > W2 > ■ ■ ■ > wt > ■ ■ ■■ 

Theorem 2 (A(t F )) assumes its upper bound if for all 
k from 1 to m 

span j=1 _ d(k) \^ r(k ^(t F )) =E(a k ), (9) 

and its lower bound if for all k from 1 to m 

span J=1 ... d{k) \^ r{k ^(t F )) =E(a N - k+1 ), (10) 

where d{k) = diva E{ak), i.e., the dimension of the 
eigenspace belonging to the eigenvalue au, and r(k,j) = 
d(l)+--- + d(k-l)+j. 

Thus, to attain the kinematical maximum for an observ- 
able whose eigenvalues are distinct with multiplicity one, 
we need to find a unitary transformation that simultane- 
ously maps the initial pure state with the largest proba- 
bility w\ onto the eigenspace corresponding to the largest 
eigenvalue of A, the initial state with the second largest 
probability W2 onto the eigenspace corresponding to the 
second largest eigenvalue of A, and so forth. If A is a 
projector onto a subspace of dimension d, then attaining 
the kinematical upper bound requires finding a unitary 
transformation that maps the d initial states with the 
d largest probabilities onto the eigenspace corresponding 
to the eigenvalue one of A. 



IV. DYNAMICAL REALIZ ABILITY OF THE 
KINEMATICAL BOUNDS 

The kinematical bounds derived above are clearly dy- 
namically realizable if every unitary operator in U(N) is 
accessible from the identity operator / via a path that 
satisfies the dynamical law (|^), i.e., if the system is com- 
pletely mixed-state controllable. 

For control-linear systems, complete controllability can 
easily be verified numerically. Note that combining eqs. 
(§) and (|) leads to 

d M 
iK-g- t U(t, t ) = H U{t, t ) + ^ frn(t)H m U(t, t ) (11) 

m— 1 

Setting x(t) = U(t,to) and 

Xm(x(t)) =- % -H m U{t,t Q ), m = 0,---,M, (12) 
equation ([ll]) becomes 

^ = X a (x(t)) + fm(m m (x(t)), (13) 

m— 1 

which defines a control system on the Lie group U(N) 
of a type studied by Jurdjevic and Sussmann H. From 
their results, the following simple algebraic conditions for 
complete controllability of control-linear, non-dissipative 
quantum systems have been derived^]: 

Theorem 3 // the total Hamiltonian is given by 
where f m are independent bounded measurable control 
functions, and dim7Y = N < oo then a necessary and 
sufficient condition for the system to be completely con- 
trollable is that the Lie sub-algebra Lo ofu(N) (the skew- 
Hermitian matrices) generated by Ho,---,Hm has di- 
mension N 2 , or equivalently, that the ideal £q ofL(U(N)) 
generated by Hi , • • • , Hm has dimension N 2 — 1 . 

If the system is not control- linear, i.e., the Hamiltonian 
depends in a nonlinear way on the control functions f m 
then there is in general no simple algebraic condition to 
verify controllability. 

For systems that are not completely controllable, dy- 
namical realizability of a particular kinematical bound 
depends on the set of unitary transformations U(t F ,to) £ 
U(N) that are accessible from the identity / 6 U(N). 
More precisely, it depends on whether the intersection of 
the set of dynamically accessible target states and the set 
of states for which the expectation value of the chosen ob- 
servable assumes the kinematical bound is empty or not. 
Since the set of dynamically accessible target states con- 
sists of all density matrices satisfying (||), where U(t F ,to) 
is a unitary transformation accessible from the identity 
in U(N) via a path that satisfies the equation of motion 
(^) , a crucial step towards answering the question of dy- 
namical realizability is to determine the accessible sets 
for the associated control system on U (N) . 
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V. DYNAMICALLY REALIZABLE BOUNDS 
FOR DECOUPLED SYSTEMS 

Among the systems that are obviously not completely 
controllable are those comprised of non-interacting sub- 
systems. We shall refer to these systems as decoupled 
systems. 

Let us first consider a control-linear Hamiltonian sys- 
tem with a single control, 



H(f(t)) = H + f{t)V, 



(14) 



where H is the internal Hamiltonian of the unperturbed 
system and V defines the interaction with the control 
field f(t). 

In this case, the system is decoupled if there exists a 
basis B for the Hilbert space H such that Hq is diagonal 
and 



V = V x © V 2 = 



Vi 








v 2 



(15) 



Let Hi and H 2 be orthogonal subspaces of H such that 
H = Hi © H 2 and each V$ maps Hi to itself, 

H = Hi®H 2 , V t :7U-> Hi, * = 1,2. (16) 

It immediately follows that H(f(t)) is block-diagonal, 



H(f(t)) =H 1 ®H 2 



hlineO 





H 2 



(17) 



and maps Hi to itself for i = 1,2. Thus, the two sub- 
spaces Hi and H 2 do not interact. Let Bi be the restric- 
tion of the basis B to the subspace Hi, Pi be the projector 
onto the subspace Hi and let N{ denote the dimension of 
Hi. 

Given an observable A on H, we define the restricted 
observables Ai = Pi A for i = 1,2. Note that Ai is a 
Hermitian operator on the subspace Hi, i.e., 



Ai 



PiAPi : Hi 



Hi 



1,2. 



(18) 



Let An denote the eigenvalues of A\, counted with mul- 
tiplicity and ordered 



A^ } > X£> > 



t(0 



> A 



(19) 



If pi(to) is the density operator for subsystem i, whose 
matrix representation with respect to the basis Bi is 



Pi{h) = diag(wf } , . . . , ) 



.,(<)• 



(20) 



with w± > w 2 > • • • > , we can apply theorem [l] 
to obtain bounds for the expectation value of A{ : 



Notice that the total probability for each subspace is 
less or equal to one, and that the sum of the subspace 
probabilities must equal one, i.e., 



Nj 



Pi 



= X>£ 3 <1, Pl+P2 = l- 



(22) 



n=l 



If the probability for subspace i is one, then the initial 
ensemble is restricted to this subspace and since the sub- 
spaces do not interact, the ensemble will remain in this 
subspace forever, i.e., pi = 1 for all times. In this case, 
(A(t)) = (A l (t)). 

If both subspaces are initially occupied, i.e., both pi 
and pi are non-zero, then the density operator for the 
entire space H is the direct sum of the subspace density 
operators pi(to) an d P2(to), he., 



p(to) = /5i(*o)® A«(*o) 



Pi(*o) 











(23) 



Since H maps each subspace to itself, we can conclude 



Pl(t) 








h{t) 



pit) = Pl(t) © p 2 {t) = 

for t > t and thus it easily follows that 

(A(t)) = Tr (Aipi(tj) + Tr (A 2 p 2 (t) 



(24) 



(25) 



Since (k(t) and Ai (i — 1, 2) are operators on Hi, we can 
apply (|2l|). Thus we have 

Theorem 4 Consider a decoupled quantum system as 
defined above. If po is given by fl29j), then the expectation 
value of an observable A is bounded by 



No 



(A(t)) > E^^-n+l + E-i 2 ^ 



(2) 

N 2 -n+ 



71 = 1 



n=l 



N 2 



i (26) 



(27) 



where A„ are i/ie eigenvalues of the subspace observable 
Ai, counted with multiplicity and ordered in a decreasing 
sequence. Furthermore, the upper bound is attained at 
t = tp if for k = 1, • • • , mi 



^. =li ... jd(fc) |*W )j) (* F )) = J E;(a^0: 

where d{k) — dim£ , (a{ £ 1 ' ) ) and r(k, j) — d(l) + - 
1) + J, and /or ^ = 1, • • • , m2 



spori i=lj ... )d ^ ) |*^ J . ) (*i?)) = E(a\ 



(28) 
-d(fe- 

(29) 



71=1 



at, w/iere d(f) = dimS(af } ) and r(^, j) = d(l) H hd(£- 

< (Ai(t)) < 'S^w^X^ (21) 1) +J- Similar conditions can be written down for the 
n lower bound. 
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This theorem provides improved bounds for decoupled 
systems and it is easy to see how it can be generalized 
to systems consisting of more than two non-interacting 
subsystems or control-linear systems with multiple con- 
trols. The improved bounds are dynamically realizable 
if all the subsystems arc simultaneously completely con- 
trollable. While the previously mentioned condition for 
complete controllability can be applied to each subsys- 
tem, it is not clear whether complete controllability of 
all subsystems always implies complete controllability of 
the system as a whole. However, our computations for 
several decoupled systems suggest that this is the case 
for the systems we studied M . 



VI. CONTROL COMPUTATIONS FOR A 
COUPLED PARTIALLY CONTROLLABLE 
SYSTEM 

For partially controllable systems whose dynamics can 
not be decomposed into independent subspace dynamics, 
the question of dynamical realizability of the kinemati- 
cal bounds for a given observable can not be answered 
in general. Rather it depends on the choice of the ob- 
servable. Among the many computations we have done, 
we studied a four-level harmonic oscillator model with 
unusual interaction terms: 



H 



0.5 

1.5 

2.5 

3.5 



Hi = f(t) 



10 
10 10 
10 1 
10 



The main difference of this model compared to the stan- 
dard harmonic oscillator model is that we set all the tran- 
sition probabilities equal to 1 instead of ^fn. Clearly, 
the system is not decoupled. Yet, unlike the standard 
harmonic oscillator, this system is not completely con- 
trollable. In fact, it can be shown that the dimension of 
the associated Lie algebra drops from 16 to 11 precisely 



when all the transition probabilities are equal. If only 
one of these values is changed, complete controllability is 
recovered jjj. 

Nevertheless, our control computations maximizing (a) 
the energy of the system, i.e., A = H and (b) the tran- 
sition dipole moment A = V, assuming 



Po = 



0.4 

0.3 

0.2 

0.1 



indicate that it still seems to be possible to control the 
system rather effectively even if all the transition proba- 
bilities are the same. The final yield after 20 iterations 
was 97.75 % of the kinematical maximum in case (a) and 
95.69% in case (b). These results suggest that further 
study of the dynamical realizability of the kinematical 
bounds for partially controllable systems is necessary. 
VII. CONCLUSION 

We have shown how the kinematical constraint of uni- 
tary evolution for non-dissipative quantum systems gives 
rise to kinematical bounds on the optimization of arbi- 
trary observables and established general criteria for the 
dynamical realizability of these kinematical bounds. It 
has been shown in particular that the kinematical bounds 
are always dynamically realizable for completely control- 
lable systems and that improved bounds can be derived 
for decoupled systems. Finally, we have demonstrated in 
the last section that certain modifications of a completely 
controllable system may lead to a loss of complete con- 
trollability, but that even in such a case the kinematical 
bounds may still be approximately dynamically attain- 
able. In latter case, further investigation of the structure 
of the Lie algebra associated with the given control sys- 
tem is necessary to determine whether a particular kine- 
matical bound for a partially controllable system can be 
dynamically attained. 
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